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Abstract 



j^ I A simple model is considered for an open system consisting of an aggregation of 

>-C ■ magnetic particles (like greigite) in the presence of a magnetic field (H), and inter- 

im . acting linearly with a bath of 3D harmonic oscillators. Using the Feynman- Vernon 

formalism, as given in Weiss (termed FVW), the time-evolved reduced density ma- 
trix ( after eliminating the bath d.o.f.'s), is examined for environmental decoherence 
j^ ■ as defined in the FVW formalism. While decoherence is usually positive for most 

two-way couplings with the enviroment, it is found that a three — way interaction 
_ involving the system plus bath plus H-field all together, can facilitate a reversal 

^ ■ of sign of this quantity ! This may have implications for quantum coherence based 

O . phenomena on the origins of life. 
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1 Introduction : Preliminaries 

o; 

O ■ Decoherence is a ubiquitous phenomenon [1] which arises when one attempts to extend 

quantum theory to macro systems (and even to meso systems ). Since environment has 
a crucial role for a quantum system, the relevant theoretical framework for the study of 

^ ■ decoherence is a theory of open quantum systems which treats the effect of an uncon- 

trollable environment on the quantum evolution. The concept was originally developed 
to incorporate the effect of friction and thermalization in a quantum formalism, but a 
correct perspective of decoherence effects requires its time scale to be much shorter than 
typical 'relaxation' phenomena. Rather ironically, decoherence is not so much about how 
a bath d.o.f. affects the system, as viceversa, the latter revealing information on the state 
of the system. The relevance of decoherence is especially acute for quantum information 
processing tasks where the coherence of many quantum systems must be maintained for 
a long time. To explain decoherence in a nut-shell, one may crudely assume that the 
interaction between the quantum system and environment is short-ranged enough to be 
described by scattering theory for asymptotically 'free' states before and after the inter- 
action, so that the detailed anatomy of the collision dynamics itself does not come into 
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the picture. Let p be density operator for the system which, for simphcity, interacts with 
a single environment d.o.f. represented hy \il) > . Then it can be shown that the diagonal 
elements of p remain unchange by interaction, i.e., p^^ = pmm ■ But the off-diagonal 
elements p^„ after interaction are no longer equal to pmn before interaction. Rather, they 
get suppressed by a factor < il)'^\%l)'^ > whose magnitude is always less than unity. And 
since the off-diagonal elements are the repository of quantum effects, this suppression is 
tantamount to loss of ability of the system to show quantum behaviour due to the inter- 
action with the environment's quantum d.o.f. [It arises from the correlation between the 
system plus environment particle]. So, if the environmental d.o.f. is integrated out, this 
ability is lost. Alternatively, the interaction is an information transfer from the system 
to the environment. The more the overlap < ip^lip"^ > deviates from unity, the more an 
observer could in principle learn about the system by measuring the environmental par- 
ticle ! But this measurement is usually never made, so this alternative principle explains 
that the wave-like interference phenomenon characterized by the coherence vanishes as 
more information discriminating the distinct " particle-like" system is revealed. 

Against this background, we may regard that decoherence that is inherent in the 
decay of an open quantum system into the surroundings, leads to a destruction of phase 
correlations within the system. The issue of decoherence of such a system due to its 
interaction with a thermal bath of harmonic oscillators was first formulated by Feynman 
and Vernon [2] and addressed subsequently by many workers, notably Grabert-Schramm- 
Ingold [3] as well as Caldiera-Leggett [4]. The result of many of these studies has been 
incorporated in a comprehensive book by Weiss [5] which will be taken as a standard 
reference for the present investigation. We shall freely use the results of ref [5], often under 
the name FVW, and in the same notation where possible without explanation. Specifically 
we consider an 'open' system consisting of an aggregation of magnetic particles (such as 
greigite) in the presence of a magnetic field (H), and interacting linearly with a bath of 3D 
harmonic oscillators. While several types of two-way interaction of the system with the 
bath has been adequately covered in FVW, the possibility of new forms of interactions 
when an external magnetic field is present, and their effect on the resulting decoherence, 
is the subject of this paper. 

In Section 2 we collect the relevant formulae from ref. [5] under the further assumption 
of slow time variation, so as to simplify the essential formulae, as well as extend them 
from 1 to 3 dimensions before identifying the decoherence rate. In Section 3 we introduce 
the H-field and examine the effect of a new form of interaction-a three — way coupling 
involving the system, bath and field. To the best of our knowledge such a direct 3-way 
coupling has not been considered in the literature [5], a novel feature being that the 
calculated decoherence has a negative sign; its origin is of course accounted for, and its 
order of magnitude estimated. Section 4 discusses the significance of this result and its 
possible interpretation on the lines of a Froehlich-like mechanism [6]. (In Appendix A we 
re-derive a formula for the problem 3-11 of Feynman-Hibbs [7], since a term appears to 
be missing in it). 



2 Essential Results for System plus Bath Coupling 

From Eq (3.11) of Weiss[5], the System-Bath Hamiltonian has the form 

V{q, x) = V{q) + \i:^m^ujl{x^ - ^f (2.2) 

This Hamiltonian serves two purposes : A) As the Euchdean Lagrangian (with imaginary 
Matsubara time [5])in quantum statistical mechanics for the calculation of the density 
matrix under thermal equilibrium conditions a la Feynman-Hibbs [7] and Weiss [5]; B) as 
the Hamiltonian corresponding to a real-time Lagrangian (obtained by changing the sign 
for the potential term V) which participates in the time-evolution of the density matrix 
for the joint system-bath complex under non-equilibrium conditions, using real time path 
integral techniques [4, 5, 7]. The Matsubara time r is related to the real time t by t = — ir, 
and varies in the range 

< r < ^/3 = h/kT (2.3) 

The global density matrix H^(0) for the initial state t=0 is assumed for simplicity to be in 
a factorized form in which the system and bath are decoupled [4,5], and the unperturbed 
bath is in thermal equilibrium, given by (see eq (5.10) of [5]): 

1^(0) = p(0) (g) Wr{Q)- WR{<d) = Z^^ exp{-/3i/fi} (2.4) 

Here p is the reduced density matrix for the system whose time evolution p{t) is the 
primary object of interest. To that end note that the global density matrix W{t) evolves 
from a state W^(0) at t=0 as (see eq.(5.1) of [5]) 

Wit) = exp{-iHt/n}W{0) exp{+iHt/n} (2.5) 

We now skip most of the formalism which is outlined in sections (5.1-2) of Weiss [5], 
and focus on the time evolution p{t) of the reduced density matrix after eliminating the 
bath d.o.f.'s. This quantity can be expressed in terms of the Feynman- Vernon influence 
functional Fpv as follows [5] 

PiQfq'pt) = jdqdq[j DqDq'eM^Ssiq] - Ss[q'])/n]FFv[q.q'] (2.6) 

Fpv in turn is made up as a product of two functionals F, F* : 

FFv[q,q] = / dxfdxidx[Wji[xi, x'^F[q; XfXi\F*[q'; Xfx'^ (2.7) 

where 

F[q;XfX,] = f Dxi.)exp[t{SR[x]) + Sj[x,q])/h] (2.8) 



and a similar expression for the (complex conjugated) quantity F*[q'; Xfx'^. The quantity 
WR[xi,x'i\ in eq(2.7) which also appears in eq (2.4), is the canonocal density matrix for 
the bath in thermal equilibrium, and its path integral over the (imaginary) Matsubara 
time r is (see eq (5.17) of [5]) a product of modes Wr^ where (see eq (5.17) of [5]) 

WRa[Xi,x[] = Caexp{ . " " [(x|, + xfj COsh (3huJa - 2XiaX[J} (2.9) 

Zii smn ptiUoi 



Similarly the quantity F[q\ XfXi\ is expressible as a product of its modes F^ whose path 
integral over real time [7,5] works out as [5] 

1 

Fa[Xfa,Xia\ = Ca{t) exp{-(t)a[q; Xfa, Xia]} (2.10) 

where Ca{t) is a constant, and 0q, is given by eq (5.19) of ref [5]. However we prefer to give 
a simplified expression for 0q, and a corresponding quantity 0'„ arising from F*[q'; Xfx'j] of 
eq (2.7), based on an (approximate) time independence of g and q', which may be adapted 
from eq.(A.7) of Appendix A. The result for (pa, mostly in the notation of ref [5], is 
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+ -i^|(-!i!i^i (2^11) 

The corresponding quantity 0^ which comes with a negative sign, arising as it does from 
F*[q';Xfx'j], is given by a similar expression to (2.11), except for 

Xia => x\^\ q^q' 

The rest of the calculation which involves the integration over the bath d.o.f.'s, is 
straightforward though lengthy. Following Weiss [5], introduce the variables 

V '^[^iat Vial -^jq, it Xj^, yA.LZj 

The integration in eq(2.6) for the reduced density matrix p, over the variables z,y,Xf, all 
standard gaussian variables, is best carried out successively in this very order. The result 
for the related Feynman- Vernon influence functional Fpv[(l,<l'], eq.(2.7), which subsumes 
all the pre-exponential factors a la ref [4], is 

FFv[q,q'] = exp{-SFv[q,q']/^} 

Spvhq'] = W^^^^^COth{PhUa/2){l-COSC0at) 

- lv^^(g^-f)(,_g^!^) (2.13) 

2 V n,rnaUji Ua 

From (2.13) we can identify the decoherence factor [5] as the flrst term S^^^ of SFv[q,q'] 
in (2.13), which is seen to be positive.lt can be recast as an integral over the "spectral 
density " deflned by [4] 

2 ^ rriryUa 



which allows it to be expressible as [5] 

1 r°° 
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5(^) = - / J(cu)^coth(/3hto/2)(l - cosut) (2.15) 



where q^ = q — q' is the spatial separation of two different locahzed states [4]. The 
decoherence rate 'jdecoh is just the time derivative of S^^^ /h, which simphfies in the regime 
of high temperature (/3 = l/ksT) for the ohmic case J{uj) = rju as [5] 

Wo /+°° , u sin ut 
Idecoh = -^:te / "^^ — 2 (2-16) 

where we have inserted an infrared cut-off frequency wq which, ahhough not needed here, 
will be found useful for later purposes (see below). We also note that the damping rate 
is negligible compared with the decoherence rate [5]. The integral is now evaluated as a 
contour integral to give 

2 

(IT) 

Idecoh = -^exp{-tcjo} (2.17) 

ph 

We now turn to a generalization of this ID formula to 3D, so as to include the effect of 
an external magnetic field. 

3 Effect of H-field on Decoherence 

The generalization of the foregoing result to 3D is merely a matter of replacements like 

qq' ^ q- q'; qXia ^ q ■ Xj^ 

in the Hamiltonian (2.1) and the potential term (2.2). More interesting coupling structures 
arise with the introduction of an external (magnetic) field. While pairwise interactions do 
not produce any new structures, there is now the possibility of three — way interactions 
which have more interesting features. The simplest form of the latter for a spin-half 
magnet of moment /iO" is 

V3 = — -^i/io- X H ■ c„Xi„ (3.1) 

where the factor i associated with fia meets the hermiticity requirement. And the factor 
in front ensures dimensional homogeneity with the formalism of Section 2. Note the 
correspondence q 4^ Q where 

Q = ^ij^d X H (3.2) 

The formalism now goes through just as in Section 2 (adapted to the 3D form) as indicated 
above, except for the replacement of q^ = {q — q'Y i'^ (2-15) by 

Q^ = ^{^\l^f2H^ (3.3) 

making use of the result {a x H)^ = 2H^. Note that the presence of the factorji}^ in 
(3.3) renders the quantity negative as a whole, thus reversing the sign of the decoherence 
factor! This brings out the central result of this exercise: The 3-way coupling involving 
the magnetic field gives a contribution to the decoherence with a reversed sign. As to 
the magnitude of the effect vis-a-vis the reference value (2.17), the counterpart of (2.16) 
reads as 

H ''■'^VQo /"^°° T 2ujjjUJ sin ut 

idecoh = _t2o / ^-,2 1, ,2\2 



Tih /3 J -00 (a;2 + co^ 
^\ = if^rH'/h' (3.4) 



where the effect of the infrared cut-off frequency appears (even more strongly) in the last 
factor. Evaluation of the integral as a contour integration over the double pole finally 
gives 

7.1o. = S^^exp{-tc.o} (3.5) 

The value of t in front of the exponential on the rhs is t ~ ^o^-, so that the ratio of the 
two decoherence factors (3.5) to (2.17) finally works out as 

Idecohhdecoh ~ i^^u/^l (3.6) 

The factor i^ brings out the reverse sign, while the smallness of cuq wrt uh (note that ojq 
is a small infrared cut-off, while ooh is a finite quantity) testifies to the enhancement of 
the (coherence) effect. 

4 Discussion 

This work was motivated by the need to study the role of a magnetic field in the context 
of origin of life scenarios [8]. While referring the interested reader for details of the 
precise mechanism involved therein, suffice it to say that the principal object behind the 
introduction of a magnetic field is its role in bringing about coherence in a system whether 
it is directed against the dissipating effect of temperature or for generating an environment 
conducive to a quantum scenario. Taking a simple algebraic viewpoint, the respective 
signatures of coherence / decoherence must be of opposite signs in whatever definitions 
are employed for their respective measures. Viewed in that light, the result (3.6) is fully 
consistent with the purported introduction of a magnetic field in the interaction of an 
open system with its environment, viz., to reduce the effect of decoherence by a reversed 
sign of the contribution of the magnetic field, taking advantage of a special 3-way coupling 
involving the system, environment and the (external) magnetic field. As to the relative 
magnitude of the factor Uh vis-a-vis Uq in Eq (3.6), the former, being a finite quantity, 
easily exceeds the latter which is a small infrared cut-off. Somewhat different approaches 
[9, 10] directed towards avoiding decoherence for open systems by adjusting external 
controllable parameters (such as a magnetic field) seem to to converge towards similar 
conclusions, although employing different strategies. 
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Appendix A. 

Here we outline some simplified expressions for the quantities (pa and 0'^ of Section 2, 
based on an elementary derivation of the classical action for a harmonic oscillator in the 
presence of a constant force /. To that end we first put on record the more general result 
for a time- dependent force f{t), vide problem 3-11 of Feynman-Hibbs [6], mostly in their 



notation, from which the result for / = const follows as a particular case. [An added 
reason for this exercise is the appearance of an extra term in the classical action, over and 
above the result listed in this classic book [7], in the hope of obtaining a verdict on its 
bona fides from interested experts]. The classical action for this problem is given by [7] 

S^i = / dt\-mx^ - -mwV + xf{t)] (A.l) 

Jo 2 2 

which can be transformed via the equation of motion as 

/■^ 1 1 

Sd = / dt[-xf{t)] + -m[XbXb - XaXa] (A. 2) 

Jo 2 2 

where Xa and x^ are the initial and final positions and their dots are the corresponding 
velocities. Now the solution of the equation of motion may be written as 

X = c exp iut + c* exp —iut -\ of(t) (A. 3) 

m D^ + co^ 

where the last term - the particular integral P/-is obtained through the following steps 

^ -e^'^* f e-'^'' f{t')dt' + C.C. 



2imu Jo 

= -^e'^^'F^it) —e-''^'F*{t) (A.4) 

2imuj 2imuj 

thus defining the integral function F^{t) and its complex conjugate F*{t). The constants 
c and c* are determined from Xa = c + c* and 

Xb = ce'^^ + c*e-'"^ + -^[e'^'^FJT) - e-'^^ F*iT)] 

2imuj 

Straightforward substitution in (A. 2) and simplification gives 



Sci = 7r~- ;7;[(a;a + ^fe)coswT-2a;aa;bH f{t)smujtdt 

2sinioI mu Jo 

2r r^ 2 

+ ^ f{t)smu{T-t)dt-—— (A.5) 

moo Jo 171^00-^ 

f dt f dsf{t)f {s)smco{T-t) sinus] + S{EXTRA) (A.6) 

Jo Jo 



X 

which agrees with Feynman-Hibbs [7] as well as Weiss [5], except for the term S{EXTRA) 
given by 

S( EXTRA) = ^— / dt f dsf{t)f{s) smuj{t - s) (A.7) 

2muj Jo Jo 

Unfortunately we are unable to write off the last term, eq (A.6), from the complete 
solution given in [7]. We also note that it does not figure in Weiss [4]either, although 
we could not detect any error in our (repeated) calculations. Its source can be traced to 
the integral term in (A. 2), wherein the complementary part of the x-solution, eq (A. 3), 
contributes an amount 

^[cF:iT) + c*FUT)] 

which remains after identifying all the terms of prob 3-10 in Feynman-Hibbs [7]. 
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